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Abstract 

The electric dipole moments (EDMs) of electron and nucleons are promising 
probes of the new physics. In generic high-scale supersymmetric (SUSY) scenarios 
such as models based on mixture of the anomaly and gauge mediations, gluino has 
an additional contribution to the nucleon EDMs. In this paper, we studied the effect 
of the C-P-violating gluon Weinberg operator induced by the gluino chromoelectric 
dipole moment in the high-scale SUSY scenarios, and we evaluated the nucleon and 
electron EDMs in the scenarios. We found that in the generic high-scale SUSY 
models, the nucleon EDMs may receive the sizable contribution from the Weinberg 
operator. Thus, it is important to compare the nucleon EDMs with the electron 
one in order to discriminate among the high-scale SUSY models. 


1 Introduction 


The standard model (SM) is established by the discovery of the Higgs boson with mass of 
125 GeV at the LHC run 1 [1-3]. However, there are several reasons to expect that the SM 
must be an effective theory of a certain full theory; no candidate of the dark matter (DM), 
no reason of the gauge anomaly cancellation, and so on. Supersymmetry (SUSY) is one 
of the attractive extensions of the SM. The lightest supersymmetric particle (LSP) is the 
candidate of the DM. The gauge coupling unification is improved due to the additional 
matters, that is, the SUSY partners of the SM particles. This unification may imply that 
the SM gauge groups are embedded in a larger gauge group, such as SU( 5), S'O(IO), and 

e 6 . 

However, the LHC run 1 has also reported that there is no signal of new physics 
around electroweak (EW) scale. Besides, the observed Higgs boson is too heavy in the 
minimal supersymmetric standard model (MSSM) if the SUSY particle masses are smaller 
than 0(1) TeV. It is needed to introduce the large quantum correction to the Higgs mass 
or the additional tree-level contribution. Several extensions of the MSSM are proposed; 
introduction of additional vector-like matters [4], specific mass spectra (large A-term or 
Next-to-MSSM) [5], and high-scale SUSY scenarios [6-8]. 

In the high-scale SUSY scenarios, sfcrmions have masses around 10 2 TeV, while the 
gaugino masses lie around several TeV. This mass spectrum leads to the fascinating results: 
the SUSY flavor and CP problems are eased due to heavy sfermions [9], and the neutral 
wino behaves as the LSP with mass of several TeV, which is favored in the thermal DM 
scenario [10]. Recently, the sfcrmion flavor structure is focused attention on in order 
to survey these models by using indirect searches [11-13]. In the grand unified theories 
(GLITs) based on the high-scale SUSY scenarios, the specific mass spectrum yields the 
several features: the gauge couplings unify at the GLIT scale with higher accuracy [14], 
and the dangerous proton decay via color-triplet Higgs exchange is suppressed due to the 
heavy sfermions [12,15]. 

The simplest model for the high-scale SUSY scenarios is based on the anomaly me¬ 
diation [16,17]. On the other hand, the generic models may include the gauge mediated 
contribution to the SUSY-breaking terms [18,19]. In the extensions, the pattern of gaug¬ 
ino masses differs from the simplest model. The vector-like multiplets for messengers are 
naturally introduced since they may obtain masses proportional to the gravitino mass via 
the Giudice-Masiero mechanism [20]. If we assume that the vector-like multiplets are in 
SU( 5) multiplets, the gauge coupling unification is maintained in these models. Thus, 
the extensions should be considered equally to the simplest model. 

The electric dipole moments (EDMs) are important to investigate the additional CP 
violation in the SUSY breaking terms. In the SM, the EDMs for fundamental fermions are 
small [21], and thus the EDMs have high sensitivities on the new physics. In high-scale 
SLISY breaking scenarios, the generic CP-violating phases are still allowed thanks to the 
heavy sfermions, even if the EDMs are generated at one-loop level. However, the future 
experiments for EDM searches may have sensitivities to the high-scale SUSY models. 

In the high-scale SUSY scenarios, the dominant contribution to the EDMs in the 
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MSSM comes from the Barr-Zee two-loop diagrams, especially the chargino/neutralino 
two-loop diagrams [22], When the higgsino and wino are around a few TeV, the current 
experimental upper bound on the electron EDM has already given the constraint on the 
models. The ratios of electron and nucleon EDMs are predictive so that the measurements 
of the ratios would lead to determination of the mass spectrum. 

On the other hand, there may exist additional contribution from gluino in the extended 
models mentioned above. The physical complex phase of the gluino mass may arise from a 
relative phase of the anomaly and gauge mediated contribution to it. The additional CP- 
violating source, so-called the gluino chromoelectric dipole moment (CEDM), is induced 
by the physical phase of the gluino mass and the CP-violating couplings of gluino and 
vector-like multiplets. The nucleon EDMs are affected by the additional source since 
the gluino CEDM turns into the CP-violating gluon Weinberg operator [23] below the 
gaugino threshold. In this paper, we study the effect of the gluino CEDM contribution 
to nucleon EDMs, and then we show the future prospects for the observation of EDMs 
(electron, neutron, and proton) in the high-scale SUSY models. 

When both left- and right-handed sfermions have flavor-violating soft mass terms, 
the one-loop diagrams to the EDMs are enhanced by the heavy fermion masses. The 
contributions may be sizable when the flavor violation is 0(1) [24-26]. However, they 
are quickly suppressed when the flavor violation is small. Then, we do not include the 
contribution to the EDMs from the flavor violation in this paper. 

Current status of the EDM experiments is as follows: the bounds on the electron, 
neutron, and proton EDMs are given by \d e \ < 8.7 x lCU 29 [e cm] [27], \d n \ < 2.9 x 
lCU 26 [e cm] [28], and \d p \ < 7.9 x 10~ 25 [e cm], respectively. 1 In future experiments, there 
are several proposals [30,31]: for instance, some neutron EDM measurements may achieve 
a sensitivity of \d n \ ~ lCU 28 [e cm]. In the proton EDM measurement at COSY [32] and 
BNL [33], they may achieve a sensitivity of \d p \ rsj 10 29 [e cm]. For the electron EDM, 
the final purpose of the ACME experiment is to reach a sensitivity of 3 x 10 _31 [ecm], If 
the electron and nucleon EDMs are discovered, we may discriminate models beyond the 
SM by taking correlation among them. 2 

This paper is organized as follows: in Section 2, we introduce the high-scale SUSY 
scenarios with the gauge mediation and show that the physical phase of the gaugino mass 
appears. The complex gluino mass gives rise to the gluino CEDM, and the CP-violating 
Weinberg operator [23] is induced by integrating out the gluino held in the scenarios. 
The gluino CEDM and the Weinberg operator in the scenarios are shown in Section 3. 
In the next section, we briefly introduce our method to estimate the observable EDMs, 
in particular neutron and proton EDMs. The contributions of (C)EDMs of quarks to 
the nucleon EDMs are evaluated with the QCD sum rules, while those of the Weinberg 
operator are based on the naive dimensional analysis. We give the detail of calculations in 
Appendices A and B. In Section 5, we study the effect of gluino CEDM to nucleon EDMs. 
We evaluate the electron, proton, and neutron EDMs in the high-scale SUSY scenarios in 
the last of this section. Finally, we summarize this paper in Section 6. 

1 The proton EDM is deduced from the Mercury EDM [29]. 

2 See Ref. [34] and the references in it for the recent works. 
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2 Complex Gaugino Mass 

In the SUSY breaking sector, we assume that there is no singlet superfield. Under this as¬ 
sumption, the soft parameters are given as follows: mass parameters of scalar components 
are induced by the Planck-suppressed higher-dimensional operators, and gaugino masses 
and scalar trilinear couplings are induced with one-loop suppression by the anomaly me¬ 
diation [16,17]. In particular, the anomaly mediated gaugino masses are given by 

M AMSB = (1) 

9a 

where the subscripts a = 1-3 denote the gauge groups of the SM, U(1)y, SU(2)l, and 
SU(3)c■ /3(g a ) and m 3 / 2 denote the beta function for the gauge coupling g a and the 
gravitino mass, respectively. 

We note that there also exists the additional contribution from the higgsino-Higgs loops 
to the wino and bino masses. Below the higgsino threshold, the additional contribution 
is given as [16] 


M* h = 


9i 3 


m\ h = JLl. 


IQtt 2 5 ’ 167T 2 
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where fin and uia denote the masses of higgsino and heavy Higgs bosons, respectively. 
tan/3 denotes the ratio of the vacuum expectation values (VEVs) of MSSM Higgs bosons. 

If there exist vector-like superfields, so-called messenger multiplets of the gauge medi¬ 
ation, the soft parameters differ from the simple high-scale SUSY breaking scenario. In 
order to maintain the gauge coupling unification, we assume that messenger superfields 
are in 5 + 5 or 10 + 10 representation. The mass terms of the messenger superfields arise 
from the Giudice-Masiero mechanism [20], if the Kahler potential is given as 


K. = |$| 2 + |<h| 2 + (c*$$ + h.c.) , (4) 

where and denote the messenger chiral superfields. On the other hand, the superpo¬ 
tential W (<h, <h) may have the mass terms of the messenger superfields, 


fU(4>,$) = Mcj/Mv 


(5) 


In those cases, the mass matrix of the scalar components of and is given as 



|M$ + c$m 3 / 2 1 2 c %m\/ 2 

c & m l /2 I + c$m 3 / 2 1 2 
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Figure 1: One-loop diagram contributing to gluino mass. Solid and dashed lines corre¬ 
spond to the propagators of fermionic and scalar components of chiral multiplets <F, <3>, 
respectively. 


where the term proportional to 717 , 3/2 arises from the Giudice-Masiero mechanism. For 
simplicity, the mass parameters are re-parametrized by M and F in the last form, and dp 
denotes the complex phase of F. 

In the above models, the gaugino masses are induced by the anomaly and gauge 
mediation mechanisms at one-loop order. The anomaly mediated gluino mass is given 
by [16,17]: 

M AMSB = ^63^3/2 • ( 7 ) 

63 = — 3 + IV5 + 31Vio is the coefficient for the one-loop beta function of the SU(3)c 
gauge coupling, where IV5 and N 10 denote the numbers of pairs of 5 + 5 and 10 + 10 
representations, respectively. We choose M amsb as real for simplicity. M gmsb is induced 
by a diagram in Fig. 1 [35]: 


Ml 


GMSB 


jl L 
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F 
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( 8 ) 


where x = \F/M 2 \ and 773 ($) is the sum of Dynkin indices of the pair of chiral multiplets, 
$ and $. The loop function g(x ) is given by 

_ (1 + x) ln(l + x) + (1 - x) ln(l - x ) (c ^ 

9\ x ) — 2 ■ 

x z 

For later use, we define the real gluino mass parameter Mg and the phase of the gluino 
mass 0 as: 


Mg e^ d = M 3 amsb + M 3 gmsb . (10) 

By the chiral rotation of the gluino field 'g a —* g ,a = e _ * 075 // 2 g a , this additional complex 
phase appears in the interaction terms between gluino and messengers. In next section, 
we show that the additional complex phase gives rise to the gluino CEDM, and then the 
CP -violating Weinberg operator also arises from the gluino CEDM operator below the 
gluino threshold scale. 
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Figure 2: One-loop diagram inducing gluino CEDM. Solid and dashed lines correspond 
to propagators of fermionic and scalar components of chiral multiplets <f>, $, respectively. 


3 Gluino CEDM and Weinberg Operator 


In our models, the gluino CEDM is generated by the messenger loop diagrams (Fig. 2) 
since there exists the non-vanishing CP-violating phase in the gluino-messenger interac¬ 
tion after the chiral rotation of the gluino. The gluino CEDM d g is given as 

C~ g cedm = ~dyf b a^ l5 G%[T a ] bc g c , (11) 

where a> w = and G“ v = d^G a v - d„G“ + gaf^Gffi. [T a } bc = - if abc and f abc is 

the structure constant for the SU(3)c- 

We estimate the relevant CP-violating terms at the gluino mass scale ( Mg ) from 
those at the messenger mass scale (M mess ) by using the renormalization group equation 
(RGE) analysis. It is useful to define the dimension-six gluino CEDM operator in order 
to estimate the RGE evolution. The gluino CEDM operator O g and its Wilson coefficient 
Cg are defined as 

0}= , d } =M,g 3 C } . (12) 

By evaluating a diagram in Fig. 2, we obtain the Wilson coefficient of O g as 

Cg(M mess ) = — 32^-2 M ~ m 2 Shl ^ + 9f ^ ^ ( ' r+ ' > + B ( r+ ^ ~ ^ m+ ’ r+ ^ r -)’ 

where m 2 ± = \M\ 2 ± \F\ are the mass eigenvalues of the mass matrix for the scalar 
components of <h and <h, and r± = \M\ 2 /m\. 6 and dp are respectively the phases of the 
complex gluino mass and the off-diagonal element F of the mass matrix m|, as defined 
in the previous section. The loop functions A(r) and B(r ) are given as 


A(r) = 


2(1 


3 - r + 


2 lnr 


B(r ) = 


2(1 


1 + r + 


2 r In r 


(14) 


Now, we estimate the gluino CEDM at the gluino mass scale by using the RGEs 
between the messenger and the gluino mass scales. The RGE for the Wilson coefficient 
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Figure 3: One-loop diagram inducing CP-violating Weinberg operator. The blob denotes 
the gluino CEDM operator. 


Cg(n) at the leading order is given as 


d 

d In /i 


CM 



I 67 r 2 




(15) 


where 7 = 12Nc and Nq{= 3) is the number of colors. This anomalous dimension is 
found by substituting the Casimir invariant Cp = Nq for Cp = 4/3 in the anomalous 
dimension of the dipole operator for b —)• sg [36,37]. The gluino CEDM at the gluino 
mass scale is obtained as follows: 




Here, 63 = — 7|sm + 2| g i U i n o where the subscripts “SM” and “gluino” indicate the contri¬ 
bution from the SM particles and gluino, respectively. The exponents —3Nc/b 3 and 1/2 
are due to the one-loop renormalization-group evolution of the gluino mass Mg and the 
strong gauge coupling g 3 , respectively. 

The CP-violating Weinberg operator is induced through the gluino one-loop diagram 
in Fig. 3. The effective Lagrangian for the Weinberg operator is defined as [23] 

Cw = C w O w , O w = — i g 3 f ahc e i lu P rJ C/ A G/ A G c pa . (17) 

Here, e livprJ is the totally-antisymmetric tensor with e 0123 = +1. At the gluino mass scale, 
we obtain the Wilson coefficient of the Weinberg operator by matching the amplitude as 
follows 3 : 

Cw{M,) = ^C 5 (M 5 ) , (18) 

3 In Refs [38-41], authors considered the similar effects arising from the CP-violating coupling of 
heavy quarks. 
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Notice that quark (C)EDMs are induced as a three-loop contribution after integrating 
out the messenger multiplets. The Weinberg operator could be also induced directly after 
squarks and messenger multiplets are integrated out, not via the gluino CEDM. These 
contributions are suppressed by the masses of squarks and messenger multiplets. Thus, 
the gluino CEDM contribution is dominant for the Weinberg operator as far as gluino is 
lighter than them. 


4 Nucleon Electric Dipole Moments 


In previous section, we show the expression of the CP-violating Weinberg operator in¬ 
duced by the gluino CEDM at the gluino mass scale. The quark EDMs are also induced 
through the Barr-Zee diagrams [42], which are dominated by only the chargino and neu- 
tralino loops in the MSSM based on the high-scale SUSY scenarios [22], The quark 
CEDMs are radiatively induced from the Weinberg operator so that they are subdomi¬ 
nant. 

In this section, we summarize the RGE evolutions of the CP-violating operators and 
results of the QCD sum rules and the naive dimensional analysis in order to obtain the 
nucleon EDMs at the hadron scale (n — 1 GeV) in a compressed way. 

The CP-violating operators in the QCD sector below the gluino mass scale are given 
as 


Per = e-S^G^G^ 


32tt 2 

\ d q q(F-a)^q- 


q=u,d,s 


l 

3 


+ -wf G^G ,up G c 


^ d qm(G-a)^q 

q=u,d,s 


(19) 


Here, F pu and G pu are the electromagnetic and gluon held strength tensors, respectively, 
and we define as F ■ a = F pu a pi/ and G ■ a = G^a'^T 0 '. The dual held strength tensor 
is defined as G a ' pu = \e pupa G a pcr . The hrst term of Eq. (19) is the dimension-four CP- 
violating term, so-called QCD d-terrn. Since this operator, however, does not mix with 
the other operators, we neglect the QCD d-term in the RGE analysis. The second and 
third terms of Eq. (19) correspond to the quark EDMs and CEDMs, respectively. The 
last term is the Weinberg operator [23]. The coefficient of the Weinberg operator w is 
given as w = —g^Cw, which is evaluated in the previous section. 

In our numerical evaluation of the nucleon EDMs, we include the RGE evolutions of 
these operators between the gluino mass scale and the scale of g = 1 GeV. The RGEs 
for the Wilson coefficients at the leading order are given by Ref. [37]. The detail of RGEs 
with mixing of the quark (C)EDMs and the CP-violating Weinberg operator is given in 
Appendix A. 

Next, we show the nucleon EDMs induced by the quark (C)EDMs via the QCD sum 
rules. In Ref. [43], the neutron EDM d n is related to the quark (C)EDMs by using the 
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Figure 4: Parameter dependence of dpr(ui) with fixed sfermion mass scale (Ms = 100 TeV). 
(Left): Messenger mass M dependence with dp = 0.1257T and x = 0.99. (Middle): Phase 
6p dependence with M = Ms and x = 0.99. (Right): x (= \F/M 2 \) dependence with 
dp = 0.1257T and M = Ms- 


QCD sum rules at the renormalization scale p = 1 GeV. Similarly, the proton EDM d p 
is also associated with the quark (C)EDMs at the scale of /i = 1 GeV. We obtain the 
relation between the EDMs for light nucleons and the quark (C)EDMs as follows: 

d p = —1.2 x 10~ 16 [e cm]# + 0.78 d u — 0.20 d d + e(— 0.28d u + 0.28#^ + 0.021#,), 
d n = 8.2 x 10 -1 '[e cm]# — 0.12 d u + 0.78 d d + e(— 0.30d u + 0.30# d — 0.014# s ). 

The explicit formulae and the numerical values are presented in Appendix B. The quark 
EDM contributions to the neutron EDM are consistent with the recent result with the 
lattice QCD simulation [44]. In the following numerical analyses, we use the results from 
the QCD sum rules with # = 0 4 . 

The nucleon EDMs induced by the Weinberg operator are given by [45]: 

d*H~e(10-30) MeVtn(lGeV), (N = n,p). (21) 

This is based on the naive dimensional analysis. The sign of the contribution of the 
Weinberg operator is also ambiguous. We adopt the value djy(w)/e = 20 MeV ru(l GeV) 
as the nucleon EDMs induced by Weinberg operator in the following numerical analyses. 


5 Numerical Results 

Now, we estimate the electron and nucleon EDMs in the high-scale SUSY scenarios. 

4 If we impose the Peccei-Quinn symmetry, the theta parameter 9 is induced, and the formulae for 
the nucleon EDMs are changed (the explicit expressions are given in the Appendix B), especially the 
coefficients of d q . However, in our study, the quark CEDMs are subdominant since they are induced only 
through the RGEs, and our results are almost unchanged in each case. 
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To begin with, let us consider the parameter dependence of nucleon EDMs induced 
by the gluino CEDM. In this evaluation, we assume that sfcrmions, heavy Higgs bosons 
and the gravitino are degenerate in mass Ms, and we take Ms = 100 TeV. Once we fix 
Ms, we have three parameters; M and \F\ which are the mass parameters of the scalar 
fields of messengers, and dp which is the complex phase of F. In the following numerical 
analyses, we choose M, 6p, and x (= \F/M 2 \ ) as independent parameters. 

Fig. 4 shows parameter dependence of the nucleon EDMs djy (N = p,n ) induced by 
the gluino CEDM. In the middle panel of Fig. 4, we show the Op dependence of gluino- 
induced nucleon EDMs. In this figure, we set M = 100 TeV and x = 0.99. If Of = 0 or 
7r, there is no contribution from the gluino-induced nucleon EDM, since M gmsb and the 
couplings of the gluino-messenger interaction are real. The maximal contribution is given 
by Op ~ 0.1257r. When the real and imaginary parts of the gluino mass are comparable 
to each other, the physical phase of the gluino mass 0 is maximized. Since the coefficient 
of the Weinberg operator is proportional to sin(d + Op) /Mg, the maximum contribution 
arises when 0 + Op ~ 7t/ 2 and Mg is small. 

In the right of Fig. 4, the x dependence of nucleon EDM is shown. This dependence 
is evaluated with M = 100 TeV and Op = 0.1257r. The gluino CEDM approaches to 
the maximum as x —> 1. On the other hand, one of the scalars of messengers becomes 
massless when x — 1. In order to avoid this situation, in the following calculations, we 
set x = 0.99 for simplicity. 

Finally, we show the M dependence of nucleon EDM in the left panel of Fig. 4. In this 
figure, we set x = 0.99 and 0 F = 0.1257T. The Wilson coefficient of Weinberg operator 
(Eq. (18)) behaves as follows: 


Gw oc 


1 


■ sin(d + Op) 


M S M 
Jp sin ( 0 + °f) 


(M « M s ), 
(M s < M). 


( 22 ) 


When M Ms, the gluino mass mainly comes from the gauge mediation so that the 
Weinberg operator is highly suppressed by M 2 and also the suppressed CP phase in sin(d+ 
0p). If M <C Ms, the gluino mass is dominated by the anomaly mediated contribution, 
and thus, the Weinberg operator is suppressed by MsM. On the other hand, in the region 
M ~ Ms, the nucleon EDMs are slightly enhanced since the gauge mediated gluino mass 
is comparable to M amsb and thus the CP phase of the gluino mass is maximal. 

Now, let us compare the nucleon and electron EDMs in the high-scale SUSY scenarios, 
in which the Barr-Zee diagrams and gluino CEDM contribute to the EDMs. In the 
high-scale SUSY scenarios, all MSSM scalar particles except the SM Higgs multiplet are 
assumed to be heavy. Within the MSSM, the main contributions to the electron and 
nucleon EDMs arise from the Barr-Zee type two-loop diagram contributions [22], The 
diagrams include loops of charginos and neutralinos, and the contributions to EDMs 
are suppressed by rrif / M 2 pn 5 ■ (The electron or quark mass rrif appears in the quark and 

5 This behavior of the Barr-Zee contribution is understood by using the effective theory of wino after 
integrating out higgsino [46]. 
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M s [TeV] 


Figure 5: Electron EDM in high-scale SUSY models with iV 5 = 0 (dotted) and 1 (solid). 
Red lines correspond to the current bound; \d e \ < 8.7 x 10” 29 [ecm] [27]. Dark gray 
and light gray regions are excluded by gluino search at the LHC with N$ = 0 and 1, 
respectively. 


electron EDMs due to the chiral nature.) Since the higgsino mass fi h is a model-dependent 
parameter, we perform the evaluation of EDMs in a range of /ih/Ms e [10~ 2 ,1]. The 
CP phase of hh is also model-dependent and therefore we set the CP phase in the mass 
matrices of charginos and neutralinos to make the Barr-Zee contributions maximized 6 . 
We note that we set tan (3 = 3 since tan/3 is also model-dependent parameter 7 . 

In the following numerical analyses, we assume that all MSSM scalar particles except 
the SM Higgs boson have the same mass Ms and three parameters which are associated 
with the gluino induced nucleon EDMs are set to be M = Ms , Of = 0.1257T, and x = 0.99 
in order to study the maximal gluino CEDM effects. The higgsino mass p# is estimated 
as follows: input value for the higgsino mass is given at the renormalization scale p = Ms, 
and then, we estimate the higgsino mass at p = fin (Ms) by using the one-loop RGEs for 
higgsino and gauginos. 

The Barr-Zee contributions are estimated as follows: the input parameters for the 
chargino and neutralino mass matrices are estimated at /i — ^h(Ms), the coupling 

6 In the Ref. [22], there exist two independent phases 0i and 02, which are the combinations of 
phases of gaugino masses, gaugino-higgsino-Higgs couplings, and fi h■ We set sin0i = sin0 2 = 1 in 
our numerical analyses, unless otherwise stated. However, since there remain the two ambiguities to 
determine the nucleon EDMs. One is the sign ambiguity in the contribution of the Weinberg operator to 
the nucleon EDMs. The other is the choice of the independent phases as sin0i = sin 02 = — 1, which is 
also allowed to maximize the absolute value of the Barr-Zee contribution. There is no need to care about 
these ambiguities as far as discussing absolute values for each contributions. 

7 We do not take care of whether the observed Higgs mass is realized in the parameter set in this 
paper, though we set the Higgs mass in the Barr-Zee contributions to be the observed value. If there 
exist additional matters contributing to the Higgs mass (such as Ref. [47]), the favored values for tan/I 
to realize the observed Higgs mass may differ from those in the simple high-scale SUSY models [6,7]. 
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Figure 6: Nucleon EDMs from Barr-Zee diagrams and gluino CEDM. Neutron (proton) 
EDM is shown in left (right) panel. Dotted lines describe the case of vanishing gluino 
CEDM, that is, we set Of = 0. Solid lines correspond to the case of vanishing Barr-Zee 
contribution. Again, gray regions in each figures are excluded by gluino search at the 
LHC. In each panels, number of messengers is set to be N 5 = 1. 


constants associated with the chargino and neutralino loops are also estimated at n — 
Hh(Ms), and other couplings are estimated at the EW scale. 

First, we show electron EDM d e in the high-scale SUSY scenarios (Fig. 5). We assume 
that there is no messenger superfield (jV 5 = 0; dotted lines in Fig. 5) and that there is 
messenger superfields (N 5 = 1; solid lines in Fig. 5). The dark (light) shaded regions 
are excluded by the gluino search at the LHC [48,49], that is, M~ ole <1.3 TeV. The 
red lines correspond to the current bound on the electron EDM measured by the ACME 
experiment, \d e \ < 8.7 x 10~ 28 ecm. It is found that the future experiment for the electron 
EDM may have sensitivities to the SUSY breaking scale of Mg ~ 10 3 TeV in each scenario. 

In the high-scale SUSY models with messengers, there are several differences from the 
models with no messengers. One is that the constraint on the SUSY breaking scale be¬ 
comes severe in the models with messengers. This is because that the cancellation between 
the anomaly and gauge mediated contributions reduce the gluino mass since we choose 
parameters to maximize the gluino CEDM. Another is that the Barr-Zee contributions 
are slightly suppressed since the masses of wino and bino become large in the extended 
models. 

Next we investigate the nucleon EDMs in the high-scale SUSY models. First, we 
compare the nucleon EDMs induced by only the gluino CEDM with the Barr-Zee contri¬ 
bution. In Fig. 6, we show the nucleon EDMs induced by only the gluino CEDM (solid 
lines) or the Barr-Zee contribution (dotted lines). The neutron (proton) EDM is shown 
in the left (right) panel of Fig. 6. If the Barr-Zee contributions vanish, the nucleon EDMs 
are induced by only the gluino CEDM. The nucleon EDMs induced by only the gluino 
CEDM are almost the same even if the quark (C)EDMs are induced via the RG mixing. 
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Figure 7: Ratio of d n w/d n bz with mass of messengers to be M s (0.1 Ms) in left (right) 
figure. Shades in each figures show excluded region by gluino search at the LHC, again. 
In each figures, we set number of messengers to be N$ = 1. 


On the other hand, if there is no physical phase of the gluino mass, the nucleon EDMs 
come from the Barr-Zee contributions. The Barr-Zee contributions are suppressed in the 
region of heavy higgsino and gauginos. Since there remains the ambiguity in the relative 
sign between the Barr-Zee contribution and the Weinberg operator, the nucleon EDMs 
are not determined if there exist both contributions. However, it is found that these con¬ 
tributions may be comparable in some region, especially heavy higgsino. Thus, we need 
to include the gluino CEDM contribution to determine precisely the nucleon EDMs. 

Now, let us estimate the ratio of neutron EDMs induced by the Barr-Zee diagrams and 
the Weinberg operator. In this evaluation, we fix the phases of Barr-Zee contributions as 
mentioned above, and we also estimate the quark (C)EDMs and the Weinberg operator 
at 1 GeV with the RGEs. In Fig. 7, we show the ratios of the neutron EDMs induced by 
the Weinberg operator (d nW ) and induced by the Barr-Zee diagrams (d n bz)- Since the 
quark (C)EDMs and the Weinberg operator mix with each other, it is not obvious to dis¬ 
criminate between the neutron EDM derived from quark (C)EDMs and that derived from 
the Weinberg operator at low energy. We, however, identify the neutron EDM induced by 
the Barr-Zee contribution (d n bz) with d n defined in Eq. (20), which is dominated by the 
Barr-Zee contribution since the RGE effect is negligible due to the one-loop suppression. 
Similarly, one induced by the Weinberg operator (d nW ) is identified with d n (w ) defined in 
Eq. (21). 

In the left figure of Fig. 7, we show the ratio |d n w/dnBz|, taking the messenger mass 
to be Ms . (We will discuss the right figure of Fig. 7 below.) In the large [is limit, the 
Barr-Zee contributions are suppressed since the EDMs induced by the Barr-Zee diagrams 
proportional to df ~ as mentioned above. On the other hand, in the small Ms 

region, the neutron EDM induced by the gluino CEDM is enhanced due to the smallness 
of the gluino mass. Thus, in the region of small Ms and large higgsino mass, the dominant 
contribution arises from gluino CEDM, though such a region is constrained by the gluino 
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Figure 8: Neutron EDMs from the Barr-Zee diagrams and the gluino CEDM in the light 
messenger case (M = 0.1 Ms). Dotted lines describe the case of no gluino CEDM. Solid 
lines correspond to the case of no Barr-Zee contribution. The gray regions in each figures 
are excluded by the gluino search at the LHC. The number of messengers is set to be 

#5 = 1 - 

search at the LHC. 

In the last of this section, let us consider the case of light messengers. The light 
messengers would realized if some symmetry is imposed to the messengers. As mentioned 
in the beginning of this section, the gluino CEDM is induced by the CP-violating coupling 
of gluino and messengers since the gluino mass is dominated by the anomaly mediated one. 
If the messengers have the mass lighter than M$, the suppression of the gluino CEDM 
becomes mild. Therefore, the nucleon EDMs induced by the gluino CEDM become large 
in comparison with the heavy messenger case. 

In the light messenger case, since the gaugino masses are dominated by the anomaly 
mediated contribution, the phase of the gluino mass 6 in the case is approximately zero, 
and then the gluino CEDM is proportional to sin 9p. Thus, in the following numerical 
evaluation, we set Op = 7t/2 in order to maximize the gluino CEDM. If the supersymmetric 
masses of messengers are sufficiently light, the lightest scalars of messengers may have 
masses much lighter than the EW scale. In the following analysis, the scale of messengers 
is set to be M = 0.1 Ms in order to avoid too light scalars of messengers while x is 0.99. 

We show the neutron EDM in the light messenger case in Fig. 8. The dotted and solid 
lines respectively describe the cases of no gluino CEDM and no Barr-Zee contribution. 
In this case, the neutron EDM via the gluino CEDM become larger than in the case of 
heavy messengers, as expected. On the other hand, the Barr-Zee contributions are the 
same in size as the case of heavy messengers. Therefore, the gluino CEDM gives the 
sizable contribution to the neutron EDM. If there exist both of the Barr-Zee and gluino 
CEDM contributions, the nucleon EDMs may be enhanced or suppressed. However, we 
do not evaluate the nucleon EDMs in this case since we have the sign ambiguity for the 
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Weinberg operator in the nucleon EDMs. 

As is the case for the heavy messenger, we also show the ratio of each contribution with 
light messengers (in the right figure of Fig. 7). In this case (M < Ads), the contribution 
from the Weinberg operator is enhanced clue to the light messenger mass as mentioned 
in Eq. (22). Thus, the d n w dominates the neutron EDMs in a broad region in the light 
messenger scenario. 

In Fig. 6, we have shown that the proton EDMs induced by the Barr-Zee diagrams 
and the gluino CEDM behave similar as the neutron EDMs. The behavior of the proton 
EDM is also similar to the neutron EDM in the light messenger scenario. Since the future 
experiments of the proton EDM may have sensitivities to \d p \ ~ 10 -29 e cm, the gluino 
CEDM effect via the proton EDM may be found. 

Before we conclude this paper, we mention the uncertainty of numerical calculation. 
The error of a s (Adz) = 0.1185 ± 0.0006 gives uncertainty on the results about 0( 1)% and 
then it does not appear in the numerical analyses. Thus, the uncertainties mainly come 
from the QCD sum rules and the naive dimensional analysis which we use in order to 
obtain the nucleon EDMs from the quark (C)EDMs and the Weinberg operator. Espe¬ 
cially, the uncertainty in the naive dimensional analysis should be large so that we might 
expect larger contributions from the gluino CEDM to the nucleon EDMs. Furthermore, 
in the above, we assume N 5 = 1 and Wo = 0 for simplicity. If more messenger fields are 
introduced, larger contribution from the gluino CEDM to the nucleon EDMs is expected 
if the CP phases are aligned to be constructive. 

6 Conclusion and Discussion 

In this paper, we estimated the nucleon and electron EDMs in the high-scale SUSY models. 
Even if the gaugino masses induced by the anomaly mediation are real, the additional 
contributions to gaugino masses may give rise to the physical phases of gaugino masses. 
We show the case in the high-scale SUSY models which is based on the mixture of anomaly 
and gauge mediations. In particular, the gluino CEDM is induced by the physical phase 
of the gluino mass mass and the C'P-violating couplings of gluino, and then it generates 
the C'P-violating Weinberg operator. 

We estimated the effect of the gluino CEDM in the extension of high-scale SUSY 
scenarios, and we showed the nucleon EDMs induced by the gluino CEDM and the Barr- 
Zee contributions. The dominant contribution to the nucleon EDMs within the MSSM 
comes from the Barr-Zee contribution of chargino/neutralino loops at two-loop level, while 
one-loop diagrams of the messenger particles generate the gluino CEDM. We revealed 
that the gluino CEDM may affect on the prediction of EDMs of nucleons in the high- 
scale SUSY models, especially in the cases of the light messengers or heavy higgsino in 
comparison with Ads■ We do not determine the total EDMs for nucleons since we still have 
large ambiguities in the naive dimensional analysis for the Weinberg operator, including 
the sign. If it is determined precisely, the proton and neutron EDMs may be found to 
behave differently in the extended models, and thus, it would be more important to detect 
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nucleon EDMs and electron EDM or the ratios of them. 
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Appendix 

A Renormalization Group Equations 

The (flavor-diagonal) CP-violating operators below the electroweak scale are given up to 
dimension-six operators as defined in Eq. (19). For considering the RGE evolutions of 
these operators, it is convenient to define the quark (C)EDM operators Of (i = 1, 2) as 

Of = —■m q ee q q(F ■ a)^ 5 q , Of = —■m q qg 3 (G ■ a)^ 5 q . (23) 

e q denotes the electric charge of the quark q and g 3 is the QCD coupling constant. The 
Weinberg operator is also given by Ow as in Eq. (17). Then, by using these operators, 
the effective Lagrangian is rewritten as follows: 

2 

Asp" = E E Cf Of + C w O w . (24) 

i =1 q=u,d,s 

The relation between the coefficients in the Lagrangian in Eq. (19) and the Wilson coef¬ 
ficients is given by: 


d q = m q ee q C\ , d q = rn q Cf , w = - g 3 C w ■ (25) 

The RGEs for the Wilson coefficients at the leading order are given as [37] 

anb c = c ' r ’ (26) 

where 

0 2 f 8 C f 0 0 \ 

C = (Cf, Cf, C w ), r = 8 C F 16 C F -4 Nc 0 . (27) 

ib7r V 0 2N c N c + 2 N f -b 0 J 

Here, Nc(— 3) and Nf are the number of colors and quark flavors, respectively. C F = 
(. Nq — 1)/(2 Nc) is the Casimir invariant and bo = —llNc/3 + 2Nf/3 is the coefficient of 
the one-loop beta function for g 3 . 

B QCD Sum Rules for Nucleon EDMs 

In order to get predictions for the nucleon and electron EDMs, we have to estimate 
the contribution to them from the parton-level interactions. In Ref. [43], the neutron 
EDM d n is related to the quark EDMs and CEDMs by using the QCD sum rules at the 
renormalization scale fi = 1 GeV. Similarly, the proton EDM d p is also associated with 
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the quark (C)EDMs at the renormalization scale p — 1 GeV. The nucleon EDMs are 
given by using the QCD sum rules as 


h/v — 


-c 0 m 3 N (qq) 

X 2 

a n 


©at, (N = p, n). 


(28) 


Here, c 0 = 0.234, (qq) = —m 2 / 2 /(m u + uid) = —(0.262 GeV) 3 is the quark condensate, 
and A a r relates the interpolation fields with the proton and neutron fields, rri^ denotes 
the mass of nucleon N. Q N is calculated through the operator product expansions (OPE) 
for the correlator of interpolation fields and are just the coefficients proportional to (qq). 
Without the Peccei-Quinn (PQ) mechanism [50] for the strong CP problem, we obtain 


0 

0 


p 


n 


(Ae u m u p u - e d m d p d )x0 + (4 d u - 
(4e d m d p d - e u m u p u )x0 + (4,d d - 


dd) 


( v 

+1 

( K - 2 ? 

du) 

+1 



Ac u d u 


C-ddct 


(4 e d d d - e u d 



(29) 


Here, e q , m q , and d q (d q ) denote the electric charge for quark q, the mass of quark q, and 
the (C)EDM for quark q, respectively, y, k, and £ are parameters which relate the quark 
condensates on the electromagnetic background with (qq)] 


F = e q xF^v (qq) , 

9s {qG A v T A q) p = (qq) , ( 30 ) 

2 g s ( qG A u T A q ^ = ieqiF^ (qq) , 


where (• • •) F denotes the vacuum expectation value on the electromagnetic background. 
These parameters are estimated by Refs. [51,52], and then the values are given by y = 
—5.7 GeV -2 , £ = —0.74, and k = —0.34. p u and pd are defined as 


Pu 


pd 


m u y 29 

[i + 

m d ] 29 


d u d,i 
: m d 

d d d u 
m u 


+ 


+ 



(31) 


The parameter ml which is associated with the VEV of q(G ■ a)q is estimated by Belyaev 
and Ioffe [51]: ml = 0.8 GeV 2 . m* is the reduced quark mass defined as m~ l = m ~ 1 + 
rrQ 1 + mj 1 . The two-loop evolution of A^ is given as follows: 


Ajv( 1 GeV) 


( a s ( 1 GeV) \ 
V OL s (m c ) ) 


f a s ( 1 GeV) + || 
X V oc s (m c ) + ff 
= -0.0439 GeV 3 , 


Qis(m c ) A 25 
a s (2 GeV )) 


2 _ 41 
9 64 


a 


,(m c ) + 


IOOtt 

154 


a s (2 GeV) + 


IOOtt 

154 


_ 6 _ 

25 


117 

154 


Ajv( 2 GeV) 


(32) 
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for Atv( 2 GeV) = —0.0480 GeV 3 . Thus, we obtain the relation between the nucleon EDMs 
and the quark (C)EDMs as follows: 


d p = —1.2 x 10 16 [e cm ]9 + 0.78 d u — 0.20 d, d + e(— 0.28 d u + 0.28 d c i + 0.021d s ), 
d n = 8.2 x 10 _1 '[e cm }9 — 0.12 d u + 0.78 dd + e(—0.30d u + 0.30 dd — 0.014d s ). 

Even if the PQ mechanism works, the theta parameter is induced as 

- m\ x - d q 
find — — / y —, 

2 m„ 


(33) 


(34) 


and then, we find the OPE coefficients O n as 


O p — (4 d u — dd) + ( k — —£ H — X 


1 777 ^ 

0 n = (4 dd — d u ) + ( k — — £ H 


4 e u d x 

^ddd c u d. 


Cddd'j i 


(35) 


Thus, we find the relation between the nucleon EDMs and the quark (C)EDMs under the 
PQ symmetry as follows: 


rfPQ = 0.78 du - 0.20d d + e(-1.2d u - 0.15d d ), 
= —0.20 d u + 0.78 dd + e(0.29d u + 0.59cta)- 


(36) 
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